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Abstract 

We compute in next-to-leading order (NLO) perturbative QCD the transverse energy 
carried into the central rapidity unit of hadron or nuclear collisions by the partons freed 
in the few-GeV subcollisions. The formulation is based on a rapidity window and a 
measurement function of a new type. The behaviour of the NLO results as a function 
of the minimum transverse momentum and as a function of the scale choice is studied. 
The NLO results are found to be stable relative to the leading-order ones even in the 
few-GeV domain. 
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Below the transverse momenta of observable jets in hadronic collisions, pr < 5 GeV, 
but still within the applicability domain of perturbative QCD, Aqcd "C Pt, there is a 
region of semi-hard parton production, px ~ 1..2 GeV. In ultrarelativistic heavy ion 
collisions this region is especially important in the formation of Quark Gluon Plasma at 
the future colliders BNL-RHIC and CERN-LHC/ ALICE: at high collision energies, the 
few-GeV QCD-quanta, minijets, dominate the initial transverse energy production ffl, 
|2], U at central rapidities during the first fractions of fm/c. The initial energy densities 
for further evolution of the system can thus be estimated based on perturbative QCD 
(pQCD)i. 

As work pdV is done during the expansion of the system, the initial transverse energy 
is not, however, directly measurable in A A collisions: the final state Et in the central 
rapidity unit has been estimated to be only 1/6 (1/3) of the initially produced Et at 
the LHC (RHIC) M. The final Et depends practically linearly on the initial one, so for 
making reliable estimates of the measurable Et the initial Et needs to be computed 
as accurately as possible. 

The average initial transverse energy produced perturbatively in the central rapidity 
unit AF of an AA collision at an impact parameter b can be computed as [RJ 

E* A (b,^,p , AY) = T AA (b)a(E T ) PoAY . (1) 

The standard nuclear overlap function T^yi(b) accounts for the nuclear collision geom- 
etry (Taa(0) ~ A 2 /(tcR\) for A ~ 200 [[|). The first moment of the semi-inclusive 
E T distribution, <j(E t ) P0 ,ay, is the pQCD quantity we formulate and compute in NLO 
below. The scale po is the smallest transverse momentum scale in the computation. It 
also governs the formation time of the system through tq ~ 1/po- To compute the ac- 
tual values of initial Et, Po has to be determined dynamically by introducing additional 
(nonperturbative) phenomenology, see e.g. We emphasize that in this work we do 
not discuss this but aim to show that the NLO pQCD formulation is field-theoretically 
well defined and perform a rigorous pQCD computation of cf(Et) Po ,ay ■ Therefore, 
at this level po is a fixed external parameter with no other physical significance than 
Po ^> Aqcd- The dependence of <r(i?r) on p will be explicitly studied and the rigor- 
ous NLO computation presented here then sets the stage for more phenomenological 
analyses. 

We generalize the LO formulation |J of <j(Et) to NLO by introducing a new type of 
infrared (IR) safe measurement functions || which contain both the rapidity acceptance 
and the definition of perturbative collisions. In getting from the 4 — 2e dimensional 
squared matrix elements to the physical quantities we apply the procedure by S. Ellis, 
Kunszt and Soper (EKS) §, §, ft §. 

The new results obtained in our study can be summarized as follows: A consistent 
and well-defined NLO pQCD formulation of <j{Et) P0 ,ay exists. Eventhough towards 
the few-GeV region the NLO results could grow rapidly relative to the LO, a stable 
behaviour of the NLO results is discovered in the range 1...2 GeV < po < 10 GeV. This 
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signals of the applicability of the pQCD in the semihard domain. The NLO compu- 
tation also brings in a new kinematical region not present in the LO, thus increasing 
the amount of perturbatively computable in nuclear and hadronic collisions. Pre- 
viously, as the actual NLO contributions were not known, various ad hoc if-factors to 
the LO formulation |J have been introduced in the literature. We have now rigorously 
computed the NLO contributions and analysed the obtained i^-factors. 

We define the acceptance window to coincide with the central rapidity unit, AY = 
{(y,4>) : \y\ < 0.5, < cf) < 2rr}, <fi being the azimuthal angle and y the rapidity. In a 
NLO hard scattering of partons, we may have one, two, three or zero partons in our 
rapidity acceptance. All the partons are assumed massless, so the transverse energy 
within AF is the sum of the absolute values pt% of the transverse momenta of those 
partons whose rapidities are in AY: 

E T = e{yi)VTi + e(y 2 )pT2 + e(y 3 )p T 3, (2) 
where the step function e{yi) is defined as in J|, 

/ 1 if Vi e AY 
e[Vi> ~ \ otherwise. 1 ' 

In the LO, the transverse momenta are equal in magnitude: pn = Pt2 = Pt- The 
perturbative scatterings can in this case be simply defined to be those with large enough 
transverse momentum, p T > p ^> Aq CD , or p T1 + PT2 > %Po- This generalizes to NLO 
as 

Pti + PT2 + PT3 > 2po, (4) 

where po is the same external parameter as in the LO case. We emphasize that po 
above is a fixed parameter which does not depend on AY. 

The IR safe measurement functions S2 and S3 |Q can now be written down. They are 
designed to answer the following question: what is the amount of E? in AY carried 
by the partons which are produced in scatterings where at least an amount 2po of 
transverse momentum is released? 

We require both the definition of Et and the definition of hard (perturbative) scat- 
terings to be included in the measurement functions, so for the 2^2 scatterings we 
define 

S 2 (pi,P2) = ®(pti+Pt2 > 2p )S(E T - [e{y l )p T1 + t{y 2 )p T 2]) (5) 
and for the 2 — > 3 scatterings correspondingly 

Ss(p\,P2,Pz) = Q(pti +PT2 +PT3 > 2po)5(E T - [e(y 1 )p T i + e(y 2 )pT2 + e(y 3 )pT3])- (6) 

We emphasize that these measurement functions are not mere generalizations of the 
jet cone R of the conventional jet case |], [7], to a rectangular rapidity window Ay 
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but they also carry information of the scale po which is not present in the computation 
for observable jets. One should also keep in mind that contrary to the Et of single 
inclusive jets, the initial Et we are considering here, is not a directly measurable 
quantity, especially not for nuclear collisions at RHIC and the LHC The question that 
we ask, however, is a well-defined one, and can be answered by making an IR safe 
pQCD computation analogous to that of observable jets H. 

Extending the definition of the semi-inclusive .Sp-distribution of [0] to NLO gives: 



da 



dEn 



PoAY 



da 2 ^ 2 da 2 :s 

po,AY dET 



dE T 
1 
2! 



PoAY 



1 
3! 



(7) 



in which momentum conservation correlates the transverse momenta as pn = — Pt2 
in the 2 —>■ 2 kinematics, and p^i = — (pt2 + PT3) m the 2 — > 3 kinematics. The [g£PS]2 
and [(iPS] 3 stand for dj>T2dy\dy 2 and dpT2dpT3d(fi 2 d(fi 3 dyidy 2 dy 3 , respectively. 

The divergencies present in the partonic NLO cross sections can be regulated by 
computing the squared matrix elements in 4 — 2e dimensions, resulting in an explicit 
e^ 1 and e~ 2 behaviour of the divergent terms. This was done first by R.K. Ellis and 



Sexton Using the 4— 2e dimensional squared matrix elements EKS have formulated 
the calculation of observable cross sections for jet production |J [| [7|, ||. We will apply 
their approach, based on the subtraction method, here. 

As explained in ||, cancellation of the divergencies takes place only if the three- 
parton measurement function 5*3 reduces to the two-parton one S 2 in soft and collinear 
limits. Our measurement functions in Eqs. (||) and (|^) clearly fulfil these criteria. It is 
also worth mentioning that while the transverse energy is a good quantity to compute, 
e.g. the number of gluons in AY would not make an IR safe measurement function 
without specifying when two nearly collinear gluons are to be counted as one. 

From Eq. (0) one obtains the first moment of the semi-inclusive Et distribution: 



ry/s du 

v(Et) Po ,ay = / dE T E T --=- 
Jo dE T 



PoAY 

where, after integrating the delta functions away in Eq. (|?p 



^ET)lrly + a{E T )l7ly, 



a{E T ) 2 ^ 2 



PoAY — 2] 



da 2 ^ 2 ~ 

[dFSh JdPS[2 S2iPUP2) 



and 



1 r da 2 ^ 3 ~ 
*(Et)^Iy =y j ldPS] 3 —^( Pl ,p 2 ,p 3 ). 



(8) 



(9) 



(10) 



3 



The measurement functions for the first E^-moment above are denoted by 



#j(pi,P2) 



+£(2/2) 



and 



^(Pl, P2,P3) 



e(j/i)pri + AV2)VT2 + e(l/3)PT3 Q(pn + PT2 + PT3 > 2p 



(12) 



where pri = |PT2 + PtsI- Naturally, also S 2 and S 3 fulfil the IR criteria, which ensures 
that <j(Et)p 0) ay is a well-defined IR safe quantity to compute. From the general form 
of Eqs. @ and fllCf ) we notice that, by replacing the measurement functions S2 and S3 
of Ref. by S 2 and S3 above, we can follow the formulation of the problem as given 
in Ref. [0] . As all the relevant formulae with a complete discussion of the cancellation 
of the several ~ e^ 1 and ~ e~ 2 singularities can be found in all detail in 0, we now 
proceed directly to the numerical evaluation. 

Our S 2 and S3 are azimuthally symmetric, the 2 -integrals giving a trivial factor 2tt, 
so basically only three- and six-dimensional integrals remain to be done in c(-&r)j^Ay 
and ct(-Et)c^aK) correspondingly. For 2 — > 2 kinematics, the integration limits imposed 
by the acceptance cuts in S2 can be solved analytically, allowing for a quick and accurate 
computation using a NAG library [jTTJ subroutine based on an adaptive subdivision 
strategy. 

Contrary to the 2^2 case, the kinematical cuts for the three particle phase space 
are very complicated due to the additional cuts introduced in the subtraction terms. To 
keep the counterparts of the subtraction terms as they are given in || (and which have 
2^2 kinematics), the kinematical cuts for the subtraction terms in the 2 — ^ 3 integrals 
must be strictly imposed through the measurement functions S3. These integrals are 
evaluated using a Monte Carlo subroutine of NAG. 

Since the subtraction procedure adopted here is exactly that of Refs. |5], ^ 0] for 
production of inclusive high-p T jets, we are also able to make use of certain subroutines 



in the JET program [0 of EKS. In particular, we have used the subroutines of JET 
for the functions in terms of which the squared matrix elements were expressed in 0. 
In addition, we have adopted EKS's book-keeping method of the parton flavours, and 
explicitly made use of their permutation tables for different subprocesses. 

For the NLO MS parton distributions, we use the GRV94-HO |13| as implemented 



in PDFLIB [[L4J]. The renormalization scale and factorization scale are chosen equal, 
= [i-F = A 4 - Our choice for the 2^3 terms is a generalization of the 2 — > 2 case 
/1 = x pr, 

H = x (p T1 + p T2 + Pts)/2 (13) 

where is a number of the order of unity. The scale /x is therefore IR safe in the 
same sense as the measurement functions required by the cancellation of the 

divergencies. 
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Figure 1: The first Et moment (t{Et) P q,ay m the central rapidity unit \y\ < 0.5 of the 
Et distribution (||). The results are shown for pp collisions at yj~s = 5.5 TeV (LHC), 200 
GeV (RHIC) and 20 GeV (SPS) as functions of the parameter pq. The solid curves represent 
the full NLO calculation, the dotted ones show the LO results computed with GRV94-LO 
parton distributions and one- loop a s , and the dashed curves stand for the LO computation 
with GRV94-HO parton distributions and two- loop a s . All the scales are chosen (Eq. [13]) as 
fJ- = (PTl +PT2 +pT3)/2. 



In Fig. [I] we plot the first Et moment <j{Ej) P0 ^y m PP collisions (i.e. no nuclear 
shadowing is included) in the central rapidity unit \y\ < 0.5, as a function of the 
minimum transverse momentum po for = 5.5 TeV, 200 GeV and 20 GeV. Based 
on Ref. 0, for AA collisions (A ~ 200) we expect the relevant range of po to be 
Po ~ 2(1) GeV for the LHC (RHIC), so we show p = 1...10 GeV in the figure to 
see whether the NLO results become very unstable relative to LO at the few-GeV 
transverse momenta. The difference between the two LO results at ^/s = 200 GeV is 
mainly due to the difference between the one- and two- loop a s (pi) used. At yfs = 5500 
GeV the differencies of the NLO and LO parton distributions also become important. 
The .fT-factors, defined as K =(full NLO)/LO, with LO in the denominator being 
either the dashed or dotted curve, can now be directly read off from the figure: At 
RHIC, as po is decreased from 10 to 1 GeV, the i^-factors range from 2.3 to 4.8 and 
from 1.9 to 2.2, the latter being the one evaluated against the truly LO result, i.e. the 
dotted line. At the LHC the corresponding values are from 3.0 to 4.2 and from 1.6 to 
1.5 when po is varied from 10 GeV down to 2 GeV. 



5 



\y\ <0.5 p =2.0 GeV 




Figure 2: Scale dependence of <t(Et) Po ,ay as a function of = 2/x/(pri + £>T2 + PT3) for 
the LHC (upper set of three curves) and RHIC (the three lowest curves) energies at fixed 
Po = 2 GeV. Labelling of the curves is the same as in Fig. [|. 



In light of the inclusive one-jet production, |F|, the ^-factors are expected to depend 
on the scale choice \i. In the best case scenario from the viewpoint of applicability of 
pQCD at pq ~ a few GeV, a region in ~ 0(1) would exist where the NLO results 
would be stable against variations of the scale choice, i.e. independent of in Eq. 
([Hp. This is studied in Fig. |2|, where the full NLO results and the two LO calculations 
are shown as a function of for a fixed value of po = 2 GeV. The values of the 
curves at = 1 are the same as in the previous figure at p = 2 GeV. At RHIC 
energies, the scale dependence of the NLO-results is very similar to the LO case. For 
the LHC, the NLO results do not depend much on the scale choice. A similar behaviour 
is, however, found also in the LO, so the stable behaviour cannot be attributed to a 
better convergence of the perturbation series obtained by including the NLO terms. 
The stable behaviour is due to the decrease of a s (fj,) and the increase of xg(x, /i 2 ) with 
larger values of N^. In the po- re gi° n considered here, these effects virtually cancel. The 
full NLO results seem therefore to be as (in)sensitive to the choice of the scale as the 
LO results are both at the LHC and RHIC energies at po ~ a few GeV. 

The following remarks are in order: First, consider the .fT-factor calculated with 
respect to the dashed curve in Fig. 1. We note that this ratio is quite stable, but 
presents a small increase towards smaller values of po- This signals that we are working 
near the borderline of perturbative QCD. On the other hand the X-factor calculated 
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against the truly LO result, shown by the dotted curve in Fig. p], does not increase, 
which implies that the NLO results are stable all the way down to po = 1...2 GeV. 
Second, the actual magnitude of the i^-factors partly follows from the fact that below 
the value = Po there is a completely new contribution to the minijet Et introduced 
by the NLO terms: the ^-distribution at Et < Po is empty in LO but in NLO 
this is not the case anymore. Therefore, more attention shoud be paid to the overall 
behaviour than to the absolute magnitude of the A'-factors. The presence of the new 
kinematical region in NLO allows us to conjecture that since this region would be 
present also in the order next to NLO, the NNLO terms will be suppressed. Third, 
based on the jet cone R dependence of the inclusive jet production 0] one might expect 
that the .fT-factors discussed above depend on the choice for Ay. However, as the LO 
contribution to a{E T ) po ^.y is almost linearly proportional to Ay, the Ay dependence 
of the .fT-factors should become much weaker than in the inclusive one-jet case. 

An attempt to compute ct(Et) Po ,ay in NLO with the semi-hard region included has 
also been recently presented in Ref . || . The formulation || is, however, a more direct 
generalization of the inclusive jet production || |6|, |7], ||: Ay is introduced and E T 
within Ay is defined as we have done above but a cut-off E in Et is introduced. 
Consequently, at Et < E , the Et distribution is empty both in LO and especially in 
NLO. The formulation || thus misses the new, calculable, perturbative contribution 
in this region. In addition, a cut-off imposed in E T causes a Ay dependence of p 
not present in our formulation. We believe that to maximally account for the Et 
production from minijets, and thus to correctly generalize the formulation of Ref. ||, 
the starting point must be to first fix the minimum overall transverse momentum 2p 
once and for all, regardless whether the partons fall into Ay or not. As discussed 
above, this results in more than just an extension of the conventional jet calculation. 

We have in this letter extended the LO formalism of Ref. [[| to NLO. Comparison 
against LO results, taken together with the fact that a new kinematical region is 
contained in NLO, shows that the NLO computation is a meaningful one, and should 
be taken into account when evaluating initial transverse energy in nuclear collisions. 

Acknowledgements: We thank K. Kajantie and V. Ruuskanen for discussions. KJE 
is grateful to S. Ellis, Z. Kunszt and D. Soper for several conversations regarding their 
NLO jet-program. We also thank A. Leonidov for discussions regarding Ref. and 
the Academy of Finland for financial support. 



7 



References 

[1] J.P. Blaizot and A.H. Mueller, Nucl. Phys. B289 (1987) 847. 

[2] K. Kajantie, P.V. Landshoff and J. Lindfors, Phys. Rev. Lett. 59 (1987) 2527. 

[3] K.J. Eskola, K. Kajantie and J. Lindfors, Nucl. Phys. B323 (1989) 37. 

[4] K.J. Eskola, K. Kajantie, P.V. Ruuskanen and K. Tuominen, Preprint JYFL-8-99, 
|hep-ph/9909456| , to appear in Nucl. Phys. B. 

[5] Z. Kunszt and D.E. Soper, Phys. Rev. D46 (1992) 192. 

[6] S.D. Ellis, Z. Kunszt and D.E. Soper, Phys. Rev. D40 (1989) 2188. 

[7] S.D. Ellis, Z. Kunszt and D.E. Soper, Phys. Rev. Lett. 64 (1990) 2121. 

[8] Z. Kunszt, ETH-TH-96/05, lectures presented at TASI 95, Boulder, June 1995; 
[hep-ph/9603235 . 

[9] A. Leonidov and D. Ostrovsky, FIAN-TD-24-98, |hep-ph/98114T7 . 
[10] R.K. Ellis and J.C. Sexton, Nucl. Phys. B269 (1986) 445. 
[11] NAG Fortran Library, Mark 18. 

[12] S.D Ellis, Z. Kunszt and D. E. Soper, JET version 3.4, 18 March 1997. 
[13] M. Gliick, E. Reya and A. Vogt, Z.Phys. C67 (1995) 433. 

[14] H. Plothow-Besch, PDFLIB Version 7.09, W5051 PDFLIB, 1997.07.02, CERN- 
PPE. 



8 



